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We present an improved measurement of the standard model CP violation parameter sin 2(j>\ 
(also known as sin 2(5) based on a sample of 45 x 10 6 BB pairs collected at the T(4S) resonance 
with the Belle detector at the KEKB asymmetric-energy e + e" collider. One neutral B meson 
is reconstructed in the J/ipK° s , ip{2S)K%, XciK%, VcKs, J/ipK*°, or J/ipKl CP-eigenstate 
decay channel and the flavor of accompanying B meson is identified from its decay products. 
From the asymmetry in the distribution of the time intervals between the two B meson decay 
points, we obtain sin20i = 0.82 ± 0.12(stat) ± 0.05(syst). The result is preliminary. 



In the Standard Model (SM), CP violation arises from an irreducible complex phase in 
the weak interaction quark mixing matrix (CKM matrix^. In particular, the SM predicts a 
CP violating asymmetry in the time-dependent rates for B° and B° decays to a common CP 
eigenstate, /cp> with negligible corrections from strong interaction^: 

Am _ t (b° -» Sep) - r(g° -> f C p) * . , . A . m 

A(t) = T(B0^f CP ) + T(B0^f CP ) = 201 Sm Amrft ' (1) 

where T(B°, B° — > fcp) is the decay rate for a B° or B° to fcp dominated by b — > ccs transition 
at a proper time t after production, £f is the CP eigenvalue of fcp, Am^ is the mass difference 
between the two B° mass eigenstates, and <f>\ is one of the three interior angles of the CKM 
unitarity triangle, defined as <f>\ = it — arg(—V* b V t d/ — V* b V cc [)- 

In this paper, we report an improved measurement of sin2^i using 45 x 10 6 BB pairs 
(42 ftr 1 ) collected with the Belle detectoH at the T(45) resonance in collisions of 8.0 GeV e 
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Table 1: The numbers of observed candidates (N ev ) 
and the estimated background (iVbkg) in the signal 
region for each fcp mode. 



Mode 


N cv 


AW 


J/ij{e + £-)K s {7T + TT-) 

J/4>(£ + t-)K s (TT°ir ) 

iP(2S){£ + l-)K°(TT + Tr-) 

4>(2S)(J/4>ir + Tr-)K%(ir + ir-) 

Xci(J/^)K° s (k + tt-) 
rj c (K + K-TT°)K^(7r + TT-) 


636 
102 
49 
57 
34 
39 
33 


31.2 
20.8 
2.4 
4.3 
2.3 
11.1 
8.9 


j/i>(£ + r)K*"(K%<K u ) 


55 


6.0 


j/^(£+r)K» L 


767 


307 



Figure 1: (Top) the beam-energy constrained mass 
distribution for all decay modes combined other than 
J/tpKl. (Bottom) the pg 118 distribution for B° -> 
J/ipKl candidates with the results of the fit. 

to 3.5 GeV e + at KEKlS, where two B mesons reside in a coherent p-wave state until one of 
them decays. The decay of one of the B mesons to a self-tagging state, /tag! *-e., a final state 
that distinguishes B° and B°, at a time itag projects the accompanying meson onto the opposite 
6-flavor at that time; this meson decays to fcp at time tcp- The CP violation manifests itself 
as an asymmetry A(At), where At is the proper time interval At = tcp — itag- At KEKB, the 
T(4S) resonance is produced with a boost of (3j = 0.425 along the electron beam direction (z 
direction). Because the B° and B° mesons are nearly at rest in the T(4S) center of mass system 
(cms), At can be determined as At ~ Az/(f3^)c, where Az is the z distance between the fcp 
and /t a g decay vertices, Az = zcp — ^tag- The average value for Az is approximately 200 /im. 

We reconstruct B° decays to the following CP eigenstates ]]: J/tpKg, if)(2S)Kg, XciKg, rj c Kg 
for £ f = -1 and J/i^Kl for = +1. We also use B° -» J/t/jK*° decays where K*° -> K%ir°. 
Here the final state is a mixture of even and odd CP, depending on the relative orbital angular 
momentum of the J ftp and K*°. We find that the final state is primarily £j ; the £y = — 1 fraction 
is 0.19 ± 0.04(stat) ± 0.04(syst^i. For reconstructed B candidates except J/ipK®, we identify 
B decays using the energy difference AE = E™ s — E^™ m and the beam-energy constrained 

mass M bc = {E^ m ) 2 — (p^) 2 , where E™* m , E™ s , and p ™ are the beam energy, the 
energy, and the momentum of the reconstructed B candidate in the cms, respectively. Figure 
[l] (top) shows the M^ c distributions for all B° candidates except for B° — > J/ipK® that are 
found in the AE signal region. Table 1 lists the numbers of observed candidates (N ev ) and the 
background (Abkg) estimated by extrapolating the rate in the non-signal AE vs Mb c region into 
the signal region. Candidate B° — > J/tpK® decays are selected by requiring electromagnetic 
calorimeter (ECL) and/or K\ and muon detector (KLM) hit patterns that are consistent with 
the presence of a shower induced by a neutral hadron. The centroid of the shower is required to 
be in a 45° cone centered on the K® direction that is inferred from two-body decay kinematics 
and the measured four-momentum of the J/ip. Figure [l] (bottom) shows the p§" s distribution, 
calculated with the B° — > J/ipK® two-body decay hypothesis. The histograms are the results 
of a fit to the signal and background distributions. There are 767 entries in total in the region 
0.20 < pf 18 < 0.40 GeV/c with KLM clusters and in 0.20 < pf 18 < 0.40 GeV/c with clusters in 
the ECL only. The fit finds a signal purity of 60%. The reconstruction and selection criteria for 
used channels are described elsewher«3. 

Leptons, charged pions, kaons, and A baryons that are not associated with a reconstructed 
CP eigenstate decay are used to identify the 6-flavor of the accompanying B meson; high mo 



^Throughout this paper, when a decay mode is quoted, the inclusion of the charge conjugation mode is implied. 




At (ps) 

Figure 2: At distributions for the events with q£f = 
+ 1 (solid points) and q£f = — 1 (open points). The 
results of the global fit with sin 2<f>i = 0.82 are shown 
as solid and dashed curves, respectively. 



Table 2: The values of sin20i for various subsamples 
(statistical errors only). 



Sample 


sin 2(j>\ 


/tag = # u (q = +1) 
/tag =B°(q= -1) 

J/ijK^lT+IT-) 

{cc)K% except J/^X§(7r + 7r) 
J/i;K* (K° s Tr ) 


0.60 ±0.19 
0.99 ±0.16 
0.67 ±0.18 
0.88 ±0.31 
1.14 ±0.23 
1.62 ± 1.10 


All 


0.82 ±0.12 



mentum leptons from b — ► c£~~P, lower momentum leptons from c — ► sl~~v, charged kaons 
and A baryons from b — > c — > s, fast pions from B° — > D^*^ (tt + , p + , af , etc.), and slow pi- 
ons from D*~ — > D°tt~. Using those tracks, two parameters, q and r, are assigned to an 
event. The first, q, has the discrete values q = ±1 that is +1(— 1) when -Btag is likely to be a 
B°(B°), and the parameter r is an event-by-event flavor-tagging dilution factor ranging from 
r = for no flavor discrimination to r = 1 for unambiguous flavor assignment. It is used 
only to sort data into six intervals of r, according to flavor purity; the wrong-tag probabili- 
ties, wi (I = 1,6), for the final fit are determined directly from data. The decay to exclusively 
reconstructed self-tagged channels are utilized to obtain wi using time-dependent B°-B° mix- 
ing oscillation: (N f - N S f)/{N f + N S f) = (1 - 2iuj) cos(Am d Af). Here N f and N SF 
are the numbers of opposite and same flavor events. The total effective tagging efficiency is 
ELi/K 1 - 2 wi) 2 = 0-270 ± 0.008(stat)+o;oog(syst), where /, is the event fraction for each r 
interval. 

The vertex position for fcp is reconstructed using leptons from J/ip decay and that for / tag 
is obtained with well reconstructed tracks not assigned to fcp- Tracks that form a K$ invariant 
mass are not used. Each vertex position is required to be consistent with the interaction point 
profile smeared in the r-(j> plane by the B meson decay length. The requirement enables us 
to determine a vertex even with a single track. The fraction of such vertices is about 10% for 
zap and 30% for z ta g. A proper-time interval resolution, i? s j g (At), consists of a convolution of 
four components: detector resolution for zqp and ztag, shift in the ztag reconstruction due to 
secondary tracks originated from charmed hadrons such as D + and D°, and smearing effect due 
to kinematic approximation in converting Az to At. We find broad outlier components in Az 
distributions due to mis-reconstruction, which are represented by Gaussian distributions. We 
simultaneously determine ten resolution parameters from data with a fit of neutral and charged 
B meson lifetimefl and obtain a At resolution of ~ 1.56 ps (rms). The width and the fraction of 
the outlier component are determined to be 36±| ps, and (0.06±g;g|)% or (3.1 ±0.4)% (multiple- 
or single-track case). 

After vertexing we find 766 events with q = +1 flavor tags and 784 events with q = — 1. 
Figure ^ shows the observed At distributions for the q£f = ±1 (solid points) and q£f = — 1 
(open points) event samples. The asymmetry between two distributions demonstrates the vio- 
lation of the CP symmetry. We determine sin2</>i from an unbinned maximum- likelihood fit to 
the observed At distributions. The probability density function (pdf) expected for the signal 
distribution is given by 

e -\At\/r B o 

V Bis (At,q,w h Z f ) = — [1 -^(1-2^)8^201 sin Am d At], (2) 

2t b o 

where we fix the B° lifetime and mass difference at their world average valueJ. Each pdf is 



convolved with the appropriate R(At) to determine the likelihood value for each event as a 
function of sin 2<pi : 



Pi = (1 - /oi) J [f sig V sig (At' , q, Wl ,C f )R sig {At - At') 

+ (1 - fsi g )V hkg (At')R hkg (At - At')] dAt' + /oiPoi(At), (3) 

where / S i g is the signal probability calculated as a function of p§° s for J/ipK® and of AE and 
Mbc for other modes. Vbk g (At) is a pdf for combinatorial background events that is modeled as a 
sum of exponential and prompt components. It is convolved with a sum of two Gaussians, Pbkg- 
Pol and f \ are the pdf and the fraction for the outlier component. The only free parameter in 
the final fit is sin 2(f>i , which is determined by maximizing the likelihood function L = FJ i Pi , 
where the product is over all events. The preliminary result of the fit is 

sin 20i = 0.82 ±0.12(stat) ±0.05(syst). 

The systematic error is dominated by uncertainties due to effects of the tails of the vertex 
distribution (±0.030). Other significant contributions come from uncertainties (a) in wi (io 02e)' 
(b) in the resolution function parameters (l^oii)' ( c ) m ^ ne J/^-K-l background fraction (1q oil)- 
The errors introduced by uncertainties in Am d and Tgo are 0.01 or less. 

A number of checks on the measurement are performed. Table 2 lists the results obtained by 
applying the same analysis to various subsamples. All values are statistically consistent with each 
other. The result is unchanged if we use the wi's determined separately for / tag = B° and B°. 
A fit to the non-CP eigenstate self-tagged modes B° — * £)M~7r + , D*~p + , and J/ipK*°(K + ir~), 
where no asymmetry is expected, yields 0.05 ± 0.04. 

We also measure CP violating asymmetry in B° — > n'K^ decays based on 45 x 10 6 BB pairs. 
Numbers of fully reconstructed events for rj(~f^)irirKg and for pjK^ are 27.7^5 and 45.4j^J, 
respectively. Flavor of the accompanying B meson is identified from its decay products. The 
decay rate has a time dependence given by 

e -\At\/r B o 

V sig (At,q,wi) = {l + ?(1 - 2wi)[S sm{Am d At) + C cos(Am d At)]j . (4) 

From the asymmetry in the At distribution, we obtain S = 0.27^;||(stat) ± 0.07(syst) and 
C = 0.12 ± 0.32(stat) ± 0.07(syst). 
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